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I. INTRODUCTION 


It is widely believed that the Universe started by big bang. When the temperature 
dropped to a few hundred MeV, the Universe experienced a phase transition from the 
quark-gluon plasma (QGP) phase to the hadron phase. If the phase transition was of 
the hrst order, the system went through the mixed phase of the quark-gluon and the 
hadron phases. In this case, the Universe would have become inhomogeneous. Such a 
non-uniform scenario of the early Universe was studied by several authors in terms of 
the primordial nucleosynthesis and the existence of the strange quark matter, which 
may be a candidate of the dark matter IH - [§] . The condition of forming non-uniform 
Universe depends largely on the property of the QCD phase transition as the critical 
temperature Tc and the surface tension a. Although QCD is the theory of quarks 
and gluons, it is difficult to solve its dynamics to get such information. 

Hence, we resort to large numerical simulations in the framework of the lattice 
QCD theory. The critical temperature and the behavior of the phase transition are 
investigated by the lattice simulations. In the pure gauge case, the phase transition 
is of the hrst order and the critical temperature is Tc = 0.26 ~ 0.28 CeV [QQ. On 
the other hand, the full QCD simulation provides Tc ~ 0.15 CeV and the order of the 
phase transition is under debate |^. The surface tension between the conhning and 
deconhning phases has been also investigated with the Monte Carlo simulation |^,^, 
in particular, on a large lattice as 36^ x 48 x 6 |^ . The recent calculations of the surface 
tension a provide the surface tension at cx/T^ = 0.0292 for Nt = 4 and cx/T^ = 0.0218 
for Nt = 6, where Nt is the value of temporal extents. Taking Tc = 0.26 CeV, we 
hnd = 80 MeV for W = 4 and = 73 MeV for Nt = 6. Here, the surface 
tension obtained in the lattice QCD is the value at the thermodynamical critical 
temperature. In general, it is hard to get the information on the metastable state from 
the lattice simulation, and therefore the upper and the lower critical temperatures 


r0| , for instance, cannot be obtained by the lattice QCD. 

On the other hand, the dual Cinzburg-Landau (DCL) theory was developed as 
the infrared effective theory of quarks and gluons In the DCL theory, color 

conhnement is provided by QCD-monopole condensation, which naturally appears by 
taking a particular gauge, the ’t Hooft abelian gauge |^. The recent lattice QCD 
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simulations indicate that the abelian part of gluons plays dominant roles at low energy, 
while the non-diagonal gluons are irrelevant in the suitably chosen abelian gauge, 
which is called as the abelian dominance [|^,Q . Hence, the non-diagonal gluons can 
be neglected in the DGL theory. The DGL theory was able to describe the linear 
potential between a static quark-antiquark pair and at the same time the spontaneous 
chiral symmetry breaking The crucial role of monopole condensation on 

conhnement and the abelian dominance were clearly observed in the recent lattice 


simulations |^- The deconhnement phase transition at hnite temperature is 

studied in terms of the effective potential using the DGL theory MM- We are then 
able to extract the surface tension from the bump structure of the effective potential 
as a function of the QGD-monopole condensate . 

The bubble formation was discussed by Kajino in the ultra-relativistic heavy-ion 
collisions, where QGP of order of the volume; V ~ (10 fm)^ might be formed 0. On 
the way of hadronization, the QGP phase may cool down to experience the mixed 
phase of the QGP and hadron, just as the case of big bang. In this case, however, 
the size of each bubble is small and hence we ought to take into account the hnite 
volume effect for the surface tension. 

In ultra-relativistic heavy-ion collisions, we face at another mixed phase system. 
Just after the heavy-ion collisions, the color-hux tube of hnite size might be formed. 
The interior of the hux tube is in the deconhnement phase while the outside is in the 
conhnement phase. In stead of the temperature ehect in the case of the ordinary mixed 
phase, the deconhnement region is produced by the color-electric held in the same 
manner as the vortex solution in superconductivity. We are also able to investigate 
this case in the DGL theory. Thus, in the latter coexistence system, the color electric 
held plays the essential role. 

We study in this paper the surface tensions of these mixed phase systems in the 
DGL theory. We consider the hat surface case, where the two phases are separated 
by one hat plane. We study also the various geometry ehect on the surface tension 
by considering various boundary shapes. 

This paper is organized as follows. In Sect. 2, we describe the properties of 
the coexistence system using the DGL theory at hnite temperature. The surface 
tensions are extracted for the hat surface case and also for the spherical case with 
various radii. In Sect. 3, we study the vortex conhgurations in the DGL theory at 
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zero temperature. We study the surface tensions for the flat surface case and for the 
cylindrical conhguration case. Sect. 4 is devoted to the summary of the present study 
and the discussion. 


II. THE SURFACE TENSION AT FINITE TEMPERATURE 


The dual Ginzburg-Landau (DGL) theory is based essentially on the abelian gauge 
hxing a la ’t Hooft in non-abelian gauge theory as QGD, which results in the appear¬ 
ance of the magnetic monopole. In the construction of the DGL theory, only two 
assumptions are made 

i) The abelian dominance. The low energy phenomena or the long range behaviors 
of the quark-gluon dynamics are dominated by the abelian part and the non-abelian 
part can be neglected. 

ii) The monopole condensation. The magnetic monopole is condensed due to the 
self-interaction, which lead to the color conhnement . 

These assumptions are strongly supported in the lattice gauge theory . 


A. DGL Effective Potential at Finite Temperature 


We consider the deconhnement phase transition at hnite temperature. At the 
quenched level, the DGL Lagrangian can be expressed in terms of the dual gauge 
helds and the monopole helds Xa, 


Cdgl = - d^B^f + Y. ■ B^)xa? - H\Xa\‘^ - 


a=l 


( 1 ) 


Here, ea, which are the root vectors of SU(3), dehne the magnetic charges of the 
monopole helds. The dual gauge helds are dehned on the dual manifold of the 
Gartan subalgebra, B^ = {B^, B^). The dual gauge coupling constant g satishes the 


Dirac condition, eg = dvr, with e being the gauge coupling constant |^. At zero 
temperature, the monopole condensate is given as |xa| = u at the mean held level for 
all a=l,2,3 due to the self-interaction term. We know that the Lagrangian is the dual 
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version of the Ginzburg-Landau theory of superconductor. The main difference from 
the Ginzburg-Landau theory is the existence of the Weyl symmetry of SU(3) group 
PU| , which stems from the SU(3)c gauge of QGD. 

We derive the effective potential at finite temperature by introducing the quadratic 
external held, JJ2a=i IXaP |0- The generating functional is written as 


Z[J] = 

= J DXaDBf^exp J d'^X^CoGL - J^\Xa 

We rewrite the monopole held as 


G = 


( 2 ) 


a=l 


Xa = (X + Xa)e*^“, 


( 3 ) 


where the vacuum expectation value y = |(0|xa|0)| is independent of the index a due 
to the Weyl symmetry. Here, the phase value is absorbed in the dual gauge held 
and turns into its longitudinal component. 

The DGL Lagrangian is written as 

C-DGL -jjz Iwr = -3A(x' - vY - 3Jf 

a=l 

- d^B^f + ^mlBl + - mlxl] 

a=l 

+2^[-2Ax(x 2 -n^) - Jx]Xa 

a=l 

+ ■ B^f{xl + 2xXa) - K^XXI + X% (4) 

a=l 

where ttib = V^gx and = ^j2X{3x‘^ — -|- J = 2\f\x. Here, source J is related 

to the mean held x, J = —2A(x^ — n^), through the condition that the linear term of 
Xa vanishes in Eq.(4). 

We obtain the partition function Z[J] by integrating over the quadratic part of 
B^ and Xa up to the one-loop level. The Legendre transform of the partition function 
gives the ehective action as 

r(x) = —i In Z[J] -|- J 3Jx^d^x 

= — J 3A(x^ — v'^Yd^x -|- i In Det{iD^^) + 3-i In Det{iD~^), (5) 
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( 6 ) 

(7) 


where factor 3 originates from the degrees of freedom of the QCD-monopole. Using 
the imaginary time formalism, the effective potential at hnite temperature is 
obtained as jl^ , 

poo 


Ueff(x; T) = 3A(x" - + 3^ / dk e In (1 - e"^7^/^) 

TT^ JO 

+ r dke In (1 - 
2 TT"* Jo 


( 8 ) 


where the temperature dependent terms are derived from the functional determinants 
in the momentum representation in Eq.(6) and Eq.(7). 

As for the parameters of the DGL theory for our investigation of the surface 
tension, we consider two cases. In one case, we take the parameter set so as to 
reproduce the inter-quark potential |^. We use this parameter set (case A) as 
the standard one to discuss all the features of deconhnement phase transition, since 
many non-perturbative phenomena have been discussed with this parameter set in the 
DGL theory [|T^,0,|T^]. This case corresponds to the dual superconductor of type-II, 
in which > ms- 

In the other case (case B), we try to readjust the DGL parameters in order to take 
into account the recent lattice QGD results. With lattice QGD simulations for the 
SU(3) gauge theory, the critical temperature has been studied in great detail. The 
continuum limit is given as 0, 


Tr 


ks = 0.625 ± 0.03 


with the expected systematic shift of 0.004 due to the inhnite volume extrapolation 
of the critical couplings. With the use of the string tension, = 0.42 GeV, we 
get a critical temperature of Tc ~ 0.26 GeV. Another constraint on the model pa¬ 
rameters can be obtained from the scalar glueball mass, because the QGD-monopole 
in the DGL theory should appear as a scalar glueball with the mass about 1-2 
GeV ll^. Very recently, Weingarten et al. worked out the enormous calcula¬ 
tion on the mass spectrum and the decay pattern of low-lying glueballs using the 
lattice QGD, and concluded that /o(1710) is identical to the lowest scalar glueball 
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with J'^ = Therefore, we identify the QCD-monopole with this scalar-glueball 
candidate /o(1710), and adopt = 1.71 GeV for the QCD-monopole mass at T = 0. 
There is some indications on the type of the dual superconductor; the QCD vacuum 
corresponds to the case between type-I and type-II p^j25[| . Since this part is not 
fully agreed upon by the specialists, we vary ms from 0.5 GeV to 2.0 GeV in order 
to study the parameter dependence of the surface tension. We take therefore, 
case A: A = 25, n = 0.126 GeV and g = 2.3, which provides m^ = 1.26 GeV and 
ruB = 0.5 GeV. The string tension is fixed to kg = 1.0 GeV/fm, which corresponds 
to = 0.444 GeV. 

Case B: m^ = 1.71 GeV, a/^ = 0.42 GeV and m^ is changed form 0.5 GeV to 2.0 
GeV. 


In order to get the feeling of the behavior of the free energy at finite temperature, 
we discuss here the results for case A. If we keep the coupling constants at finite 
temperature in the DGL Lagrangian as those hxed by the properties of the QCD 
vacuum at zero temperature, the critical temperature comes out to be too high, 
Tc = 0.49 GeV, as compared to the lattice results. Since we work with the effective 
Lagrangian, we could vary any of the coupling constants with temperature. In fact, 
the asymptotic free property suggests the reduction of the coupling strength as the 
temperature increases. We follow the work of Ichie et al. [|T3 and take temperature 
dependence for A as 

'Tc - aT 


MT) = A (: 


T, 


( 9 ) 


and fix the parameter a so as to reproduce Tc = 0.28 GeV. They have found a = 0.88 
and discussed the properties of various quantities at hnite temperature. We show in 
Fig. 1 the behavior of the effective potential in this case. At small T, the absolute 
minimum in Ves appears at hnite y, which is indicated by x. The minimum point of 
X moves towards the smaller value with temperature. Around Tiow ~ 0.26 GeV, the 
local minimum appears at y = 0. This local minimum becomes the absolute one at 
the critical temperature, Tc = 0.28 GeV. At T > Tc, the absolute minimum stays at 
X = 0, which physically corresponds to the deconhnement phase. 

We should make one more comment on the lower critical temperature Tiow It is 
dehned as the temperature at which a local minimum of Wfr appears at the trivial 
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vacuum x = 0- Tiow is exactly derived using the high temperature expansion 

1/2 


as 


Tiow = 2v 


Q^iTlow) 


y2X{Tion,) + 3g^) ■ 

We hud Tiow = 0.262 GeV for case A. We can check this formula through the agree¬ 
ment with the numerical calculations as shown in Fig.l. We get similar results for 
case B for the free energy. 


B. The Surface Tension on the Flat Interface at Finite Temperatnre 


In general, a mixed phase is realized for the case of the hrst-order phase transition. 
In the mixed phase, the system consists of two different phases, which are divided 
by an interface. The free energy densities of both phases are equal there. The free 
energy of the mixed phase is different from the sum of the two free energies by the 
amount Fg = aS, where S is the area of the interface and a is the surface tension, 
which is the energy per unit interface. 

In the DGL theory, such a mixed phase consists of the QGP and the hadron phases. 
The free energy of the system is estimated using the effective potential Kfr(xj T) as 


nx-,T) 



(Vx)' + Wff(x;T)l, 


( 11 ) 


where x{^) is the monopole condensate in Eq.(3). The coefficient 3 of (Vy)^ is the 
degree of freedom of monopoles originated from the Weyl symmetry for the SU(3) 
case. 

The surface tension on the planar interface can be estimated analytically, when 
the effective potential Kff(X) ^c) between two minima at is approximated by a sine 
curve 


l6|. We show in Fig. 2 and the approximate sine curve. 


V;in(x) = ^{^- cos , (12) 

where xh is the monopole condensate at the non-trivial local minimum of the effective 
potential, and h is the barrier height of the effective potential between x = 0 and 
X = Xh- Here, we set V'eff(0;Tc) = 14fr(XH;^c) = 0. In this case, the system holds 
translational invariance in the direction parallel to the interface. Hence, x depends 
only on ^-coordinate perpendicular to its flat interface, and the surface tension is 
reduced to the one-dimensional integration over the energy density. 
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(13) 


[Xsin] / dz 


'dxs 


dz 


+ Kin(Xs 


Here, we use ckoo because the flat surface case can be regarded as the hadron bubble 
with the infinite radius. The held equation for ;\:sin(^) is found to be the sine-Gordon 
(SG) equation, and can be solved analytically as m 


2a/6_ / 

Xsin[z) = ^-XHtaxi e> 


(14) 


with e = Here, 2e denotes the thickness of the boundary between the two 

phases. Substituting this solution into the expression for the surface tension in 
Eq.(13), 

4^3 


C^r>n — 


TT 




(16) 


Using the formula, we And = 114 MeV for case A. 

On the other hand, the surface tension on the planar interface is also calculated 
numerically from Ves directly. The profile of x{z) is given by the differential equation 
as 


d‘^X _ 
dz"^ 6 


(16) 


which is derived by the minimization condition of the free energy in Eq.(ll). The 
surface tension is expressed as 


^OO 



+ VMx) 


( 17 ) 


We solve Eq.(16) for x and And = 117 MeV for case A. This value is very close 
to the analytical result using Eq.(15). Hence, the use of the sine-Gordon (SG) kink 
solution in Eq.(14) is good for the surface tension. This result seems to be consistent 
with the recent lattice results, = 80 MeV for A7 = 4 and 73 MeV for A7 = 6 0. 

We study also the parameter dependence of the surface tension. In the super¬ 
conductivity, the Ginzburg-Landau (GL) parameter k is useful for categorizing the 
superconductor in terms of the magnetic response; the type I and type H supercon¬ 
ductors correspond to k < l/-\/2 and k > l/\/2, respectively. This GL parameter 
is also meaningful for the dual superconductor picture of the QGD vacuum given by 
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K = m^/\/2mB in the DGL theory. Some lattice simulations suggest that the QCD 
vacuum would correspond to the case near the border {k = 1/v^) between type-I 
and type-II In order to understand the behavior of the surface tension with 

respect to k, we investigate various cases with 0.61 < n < 2.4 in the case B. We 
list in Table 1 k, g, A, u, A(T = Tc) and v{T = Tc). Here, we introduce also 
the temperature dependence of u(T = Tc), in the same manner as the temperature 
dependence of A(T = T^). v{T = T^) is obtained by varying v so as to reproduce 
Tc ~ 0.26 GeV as A is fixed to the value at zero temperature. 

We show in Fig. 3 surface tensions cToo, which are derived by using A(T = Tc) 
and v{T = Tc) as a function of the GL parameter k. The surface tensions decrease 
with K. The surface tension is found to be lOSMeV < < 241MeV for this range 

of K. If the mass is small, the value of the surface tension is close to the value, 
= 80 MeV, obtained in the lattice QGD simulations 0. 

Goleman discussed the scenario of the bubble formation in the early Universe 0- 
In the thermal history of the Universe, as the temperature drops gradually across the 
critical temperature Tc in the early Universe, the hadron bubble formation in the bulk 
QGP phase begins at Tc. It is considered that the phase transition from the QGP 
phase to the hadron phase is completed above T/o^, because the latent heat turns 
into the shock wave emission from the hadron bubble in the beginning of the phase 
transition and its reheating effect suppresses the formation of new hadron bubbles. 

When T > T^, the system is in the deconfinement phase. Near T^, only extremely 
large hadron bubble can be formed, since the formation of a mixed phase necessarily 
needs the additional surface energy and the free energy of two bulk phases are almost 
equivalent For Tiow < T < Tc, the finite radius R{T) of the hadron bubble is 

determined by the balance between the volume and the surface energies of a bubble. 
There appears the energy difference between the interior and the exterior regions of 
the bubble, AV(T) = \ Wfr(0;T) — Vefi{xH',T) |, as a source to compensate the extra 
free energy density caused by the interface region. The sizable bubble configuration 
is allowed energetically 00 ■ From this consideration, we may assume that only 
the volume energy and the surface energy contribute to the free energy F{T) as. 


A'jr 

F{T) ~ -AV(T)—R{T)^ + a^inR(Tf, 


( 18 ) 


where the source term is approximated by the surface tension of the inhnite size 
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bubble Qfoo- In many studies, the surface tension of the hnite size bubble a{R) is 
replaced by that of the inhnite size bubble for simplicity, even for the free energy of 
the hnite size bubble. For a small bubble, its free energy would be largely different 
from F{T) evaluated by Eq.(16). For the realistic estimation on the free energy of 
the hadron bubble, we ought to consider the curvature effect on the surface tension. 


C. The Surface Tension in the Hadron Bubble with the Spherical Geometry 


We discuss, here, the radial dependence of the surface tension of hadron droplets 
in the bulk QGP phase below Tc. The free energy of a hadron bubble in the spherical 
geometry is given by 

'dx^ ^ 


F(r,T) = jdriTir^ 3 + V^(r,T) , 


(19) 


where y = y(r) is the solution of the equation of motion in the spherical coordinate, 

( 20 ) 


fPx _ ^ 1 llUff 

d^2 j. dr Q dx 


with the boundary conditions for the hadron bubble; y = xh at r = 0 and y = 0 at 
r —>• cx). 

Generally, in 1+d dimensional scalar theory, the soliton equation does not neces¬ 
sarily have a nontrivial static solution because of the Derrick theorem |^. The 1-1-1 
dimensional held equation, which corresponds to Eq.(16), has always the solution to 
stabilize the free energy because the scale transformation as r' = r]r for Frj{x,T) 
supplies a stable point on t], 

f / o \ 2 

roo / rt \ 

F,{X-,T) = 


dr 




= rj (Kinetic part) -I — (Potential part). 

r] 


( 21 ) 


On the other hand, if d > 2, there exist no nontrivial solutions, because the scale 
transformation of the radial coordinate does not make the free energy stable. For 
instance, in 1-1-3 dimensional case, which corresponds to Eq.(20). There is no stable 
point on 7] in the free energy, Fr^{x')T) as 
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( 22 ) 


roo f F) \ ^ 1 1 

F^{X;T) = I^ indr'r'^ 3 |^_^(/)j _ + v;^[x(r')]^ 

= - (Kinetic part) H—^ (Potential part). 

r] T}-^ 


Thus, there is no nontrivial static solution in the 1+3 dimensional scalar held theory, 
even if one resorts to any numerical methods. However, our aim is not to hnd the 
static solution of Eq.(19), but is to study the properties and dynamics of the hadron 
bubbles, which appear due to the thermal huctuation at the early Universe. In 
general, the hadron bubbles are unstable against the scale transformation, so that 
they should expand or shrink as the time goes [P9|] . 

Here, we concentrate ourselves on the hxed size hadron bubble with the spherical 
geometry, and study its dynamical evolution in terms of the surface tension. When 
the bubble radius is large enough, the solution of Eq.(20) in spherical coordinate is 
reduced to that of the flat surface case, since we can neglect the hrst term in the 
right hand side at i? —oo. We adopt the “sine-Gordon (SG) kink form” for the 
prohle function of the monopole condensate x(r). 


X{r) = Xh 


tan ^ 


(23) 


which reproduces the exact form of Eq.(16) at T^. As shown in Fig. 4, R is the bubble 
radius dehned at the half of the monopole condensate xh at which Wfr is absolute 
minimum and 2p corresponds to the thickness of the interface. We normalize the 
prohle function so as to satisfy x(r = 0) = Xh^ which is the non-trivial minimum of 
the effective potential. In the supercooling environment, we take the minimization 
condition by making the free energy in Eq.(19) minimum for a hxed radius R by 
varying p, in other words, by changing the inclination of the SG-kink curve. We show 
in Fig. 5(a) and Fig. 5(b) the prohle of the monopole condensate x(r) and the free 
energy density at T = 0.2795 GeV, as an example. As seen in Fig. 5(b), the free 
energy density is negative inside the bubble and positive at the bubble surface. 

Subsequently, we plot the free energies for various radii R in Fig. 5(c) at T = 
0.2795 GeV as the typical case. The free energy has a peak at the critical bubble 
radius Re, that is, a branch radius for a bubble to shrink or grow during the phase 
transition. In general, there appears the critical radius Rc at any temperature from 
Tiow to Tc. One hnds that Fig. 5(c) behaves in the same manner as the free energy 
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in Eq.(18), although now the surface term is calculated more accurately. 

We estimate the surface tension at Rc, because the creation of the hadron bubble 
with Rf. is relevant in the hadronization process. Once Rc is dehned at various tem¬ 
peratures, the surface tension of the hadron bubble is derived by the substitution of 
ctoo for a[R{T)] in Eq.(18), 

aW = + (24) 


where the volume energy AV (T) does not dependent on Rc at each temperature. 

For case A at 0.278 GeV <T< 0.2799 GeV, the radial dependence of the surface 
tension is shown at various temperatures in Fig. 6(a). One hnds a{Rc) gradually 
approaches a^o with Rc. On the other hand, a{Rc) increases rapidly with decreasing 
Rc below Rc ~ 30fm, although this behavior seems of inconsistent with Kajantie et 


a/.’s assertion on the surface tension. 

This behavior of the increase of a{Rc) at small Rc can be explained from the 
increase of the surface energy by the derivative term dx/dr. As shown in Fig. 6(b), 
the thickness of the bubble wall, which is about several fm, is found larger than that 
in the lattice study . The behaviors of the free energy at various temperatures are 
shown in Fig. 6(c). They are similar to the behaviors of the free energy in Eq.(18). 
However, quantitatively, there is a difference between the free energies with a^o and 
a{Rc) in the low energy region of interest for hadronization. As for case B, a[R{T)] 
behaves similarly as that of case A. 

We show the critical bubble radius Rc during the phases transition at various tem¬ 
peratures in Fig. 7. Hadron bubbles would be born with various sizes corresponding 
to each temperature, and then would grow up with the expansion of the Universe. 
The slight decrease of the temperature rapidly makes the possible bubble radius small. 
Since the bubble formation rate is expressed as P oc exp{—47rPc(^)^^(^)/3T}, with 
h being the barrier height between the false vacuum and the true one, the formation 
rate of the hadron bubble is drastically changed with temperature near Tc. 


III. THE SURFACE TENSION ON THE FLUX TUBE 
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We study here the surface tension in the color-flux tube system with considering 
the geometrical effect of the tube. In the relativistic heavy-ion collisions, the huge 
size color flux tube is produced between the target and the projectile. This scenario 
of forming huge size color flux tube has been pointed out by Ichie et al. . 

In the previous section, the coexistence of the two phases is brought by the thermal 
effect. On the other hand, the deconhnement phase is realized by the external color- 
electric held for the hux tube system. Hence, to begin with, one should consider 
the free energy of the QCD system with the external color-held. In the case of 
the hux tube, the deconhnement phase contributes to the equivalence of the free 
energies of each phase by the induced color-magnetic current in the conhnement 
phase. This conhguration is easily inferred by analogy with the coexisting system 
of the superconductor immersed in the magnetic held. The interior of the hux tube 
corresponds to the normal phase in the superconductor and its exterior corresponds 
to the superconducting phase. 

First, we rewrite the interaction terms of the monopoles Xa and the dual gauge 
helds Bfj^ in the DGL Lagrangian into the simple form similar to the well-known GL 
theory. The dual gauge helds = {B^, B^) are induced by the color current brought 
by each quark-antiquark pair in the hux tube. The quark charges Qa = {Ql, Qa) 
determine the ratio of Bf, on Bf, as 

/.i fX 


B,= 


B,= 


^ ]b,^ 


, 2 ’ 273 / 

' 1 1 


2’ 2^3, 


B,. 




for 


for 


for 


Qi — 

<32= I 

Qs = 


.2’ 2^3, 
1 1 


e, 


2’ 2^3, 
^ 1 \ 
Vs, 


(25) 


e. 


As for Bfj_, the relevant terms appearing in the DGL Lagrangian are ea ■ B^, which 
have cyclic property on the index a. 


(^1 ■ B ^ 0, 62 ■ B ^ 63 ■ B ^ 

G ■ = 0, €3 ■ Bf^ = —G ■ B^ = — 

ds • Bfj_ = 0 , ei ■ B^ = —62 ■ B^ = ——Bfj^, 

Here, e is the root vector of SU(3), 


for Qi - Qi, 
for Q 2 - Q 2 , 
for Qs-Qs- 


( 26 ) 
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t) 

It should be mentioned that two non-trivial held equations for Xa become equal, 
and there is one trivial solution |x| = u, which is determined by the energy minimum 
condition. As a result, we hnd the following relations of the components of y for each 
color charges of quarks. 


el = (1,0). (27) 


Cl — 



€2 = 


1 

' 2 ’ 


Ixil = V, 

IX2I 

= IXsl 

= Ixl, 

for 

Qi — Qi 

IX2I = V, 

IXal 

= Ixil 

= Ixl, 

for 

Q2 — Q2 

IXal = V, 

IXil 

= IX2I 

= Ixl, 

for 

Q3 — Q3 


(28) 


From this consideration, whichever the quark-antiquark pair we take, the DGL 
Lagrangian (1) should be reduced to, 

Cdgl = ~ + 2|((9^ - igB^)x\‘^ - 2A(|x|^ - v^f. (29) 

Furthermore, we rescale the helds and the coupling constants as 

2 1 

x^^X, A^2A, g ^ ^g, ^ ^ 

which reduce £_dgl into the GL type Lagrangian for the superconductor theory, 

Cdgl = - dyB^f + |(<9^ - igBf,)x\^ - A(|x|^ - vY- (31) 

From now on, we use this simple Lagrangian with the rescaled variables. 


A. The Surface Tension in the Mixed Phase with Color-Electric Field: 

Flat Surface Case 


To begin with, we consider the ordinary surface tension on the flat surface case, 
which can be regarded as the inhnite radius limit for the cylindrical flux tube. 

Let us set the interface yz-^Aane (x = 0). Then, the system essentially depends 
only on the x-coordinate. The differential equations are given by the minimum con¬ 
ditions of the free energy F; SF/Sx* = 0 and SF/SB = 0. Taking B = (0,F(a;),0) 
and V X B = D = (0, 0, D{x)), they are found to be 
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(jV + - 2Ati\ + 2A|xpx = 0 

VxD + 292|^|2B = 0. 


(32) 

(33) 

oo (normal phase), 


We impose the boundary condition as y = 0 and D = Ec at x 
and X = n and D = 0 at x ^ +cx) (superconducting phase). 

In this mixed phase, the color flux is regarded as the external held for the non- 
perturbative vacuum so that we take into account the contribution from the external 
held as — / EdD. Here, the induced color-electric held D is brought by the external 
held E. If E is independent from D, the free energy of the system with the external 
fields / is represented as f = fo + f — DEc, where / is the free energy contributed by 
the change of the monopole condensate and /o is that of the whole system without the 
external fields. On the other hand, the free energy in the normal phase at a; —a —cx) is 
/e = fo — E‘^/2 because D = Ec and — / EdD = — f D dD\E=Ec = Here, Ec 

denotes the critical color-electric held and is determined so as to All the gap between 
the free energy densities for x = 0 and x = as i?^/2 = Hence, one can rewrite 
fE = fo — An^ The surface tension on the flat interface is given by the difference 
between / and Je, 


^OO 


/ OO ^ ^ 

dx{f - /e) 

-OO 


roo 

/ dx 

'd‘^ 

-h 

dx 

J — OO 

2 

dx 




2\v^\x? + \\X\^ - E,D + 



(34) 


In order to make these equations simple, it is convenient to introduce dimensionless 
variables expressed as x = x/S, x = x/'v,B = B/{Ec5) and D = dB/dx = D/Ec. 
The GL parameter k is given as 


K, 


2 _ 


e 



A 


(35) 


where S is the penetration length of the color-electric held 6 = l/niB and ^ is the co¬ 
herent length of the monopole condensate as .^ = \/2/m^. As the definition, k < 1/^/2 
corresponds to the type-I superconductor and k > l/\/2 the type-H superconductor, 
respectively. 

The held equations written with the dimensionless variables are 


X 

B 


// 

// 




7;B^X -X + X^ 


= fB. 
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(36) 

(37) 



The boundary conditions are given byx = 0,-D = lata;—cx) and y = 1, Z) = 0 
at X ^ oo. The surface tension is simply expressed as 

a^= [ 5El + B\B' - 1)1 dx. (38) 

J L 

We show in Fig. 8 versus n with the parameter set B. moves from 
180 MeV to —310MeV for 0.61 < k < 2.4, when is hxed as = 1.71 GeV. 
At K = 1/v^, the QCD vacuum is just at the border between type-I and type-II of 
the dual superconductor, and the surface tension is zero, which is consistent with the 
ordinal superconductivity. 

It is a remarkable feature that the surface tension may be negative in the type- 
II dual superconductor. In this case, it is more stable that a thick flux tube splits 
into many flux tube units, because the larger interface of many flux tubes make the 
free energy of the whole system lower. Hence, the huge flux tube would be unstable 
against the splitting, which may provide the inhomogeneity in the QGP formation. 


B. The Surface Tension in the Flnx Tnbe with Cylindrical Geometry 


We consider here the color flux tube with cylindrical geometry, and study the hnite 
size effect on the surface tension in the flux tube. Setting the flux tube along the 
; 2 -axis, the color-electric fields D(r) = (0, 0, B^)) is given by the vector potential 

B = (0, B^{r), 0) using the cylindrical coordinates (r, 0, z). Taking the monopole held 
as X = x(r)e*”'^, the DGL equations (32) and (33) are expressed as 


f df \ dr 


2 1 / ^/2n ^ ^ ^3 

=« X— --B] x-X+X , 

2 \ K r 


did \ (\/2n 

^ ^ (Vh ~ J ^ ^ ^ 

with the dimensionless variables as f = r/6, x = xl"^-, B = B/{Ec6) and D = 
B' + 13/r = {B' + Blr) I Ec. 

Using the cylindrical coordinate, the snrface tension is written in the same manner 
as the hat interface case. 


a{R) = 


2tiR J 2 \k 


Eld 2 .,2 


- B\ x^ -2x^ + x^ + {b - if 2nrdr, (41) 



where R is the cylindrical radius of the flux tube. Here, we dehne the flux tube 
boundary as the point satisfying, x = 0.5, that is, the monopole condensate y takes 
the half value at the boundary. 

The integer n is the winding number related to the flux quantization, and corre¬ 
sponds to the vortex number in the superconductivity. We study the prohle function 
x{r) of monopole held with various number of n. In discussing the radius and the 
multiple winding number, it is meaningful only for type-I case, (t{R) > 0, where 
one large size hux tube is preferred. Otherwise, the negative surface tension leads 
the fragmentation of the huge hux tube into many thinner hux units. The numerical 
results for x(r) are shown in Fig. 9(a) and (b), using the parameter set, k = 0.61, 

= 1.71 GeV, niB = 2.0 GeV and v = 0.163 GeV, which corresponds to the type-I 
superconductor. In Fig. 9(b), we plot the numerical results on the hux tube radius R 
and the winding number n. Here, we add the htting curve, R = fy/n with f = 0.21 
fm. Thus, R increases proportionally with y/n. From n = 1 data in Fig. 9(a), the 
radius of one hux tube is estimated about f = 0.21 fm. As the physical interpreta¬ 
tion, the thick hux tube with the winding number n can be regarded as the gathering 
system of n hux tubes with the radius f and unite winding number, and therefore its 
cross section S is given as S' = tiR^ ~ The surface tension u{R) goes down with 

the radius R related to n, and asymptotically approaches to a^o in the hat surface 
case at i? —>■ cx) as shown in Fig. 10. 


IV. CONCLUSION 


In the dual Ginzburg-Landau theory, we have studied the surface tensions in two 
diherent systems. One is the mixed phase near the critical temperature and the 
other is the coexistence system of two phases with external color-electric held at zero 
temperature. 

The surface tension at hnite temperature is a relevant quantity for the argument 
of the bubble formation in the big bang scenario. In particular, the huge surface 
tension leads to strong inhomogeneity in the early Universe |I|], which ahects pri¬ 
mordial nucleosynthesis and the formation of the strange quark matter. Not only 
for the cosmology, the surface tension is important also for the QGP creation in the 
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ultra-relativistic heavy-ion collision. We have investigated the surface tension for the 
hadron bubbles appearing in the QGP phase using the SG-kink curve ansatz for the 
prohle of the monopole condensate. 

In the flat surface case at the critical temperature, we have hrst approximated the 
effective potential by a sine curve to investigate the mixed phase analytically. We 
have thus obtained the useful analytical expressions for the surface tension and the 
thickness of the boundary surface in terms of the shape parameters of the effective 
potential. We have then calculated the surface tension numerically, and have verihed 
the validity of the analytical formulae. The surface tension aoo is obtained as = 
117MeV for case A, where case A corresponds to a typical parameter set for type-II 
dual superconductor. 

We have studied the surface tension with varying the Ginzburg-Landau parameter 
K = m^/\/2mB- The value for the surface tension is lOSMeV < < 241MeV for 

0.61 < K < 2.4. The surface tension tends to be small for large k case corresponding 
to the type-II dual superconductor. These values should be compared with = 80 
MeV obtained by lattice QGD simulations. 

We have studied also the hnite size effect on the surface tension in the hadron 
bubble with various radii for Tiow < T < Tc. To this end, we have taken an ansatz for 
the prohle function x(r) of the monopole condensate in the radial coordinate. Using 
a hxed value of the bubble radius R, the thickness parameter p is determined by the 
minimization condition of the free energy. The free energy of the hadron bubble has 
a maximum at the “critical” radius i?c, which is a branch radius of a bubble to shrink 
or grow during the phase transition. We have estimated then the surface tension at 
the critical radius i?c, «[-Rc(7")], as the relevant case for the created hadron bubble 
We have found that a (Re) stay almost the value at ctoo down to the radius 
Rc ~ 30fm, a{Rc) tends to increase with decreasing Rc. 

We have studied also the surface tension in the mixed phase with the external 
color-electric helds at zero temperature. Such a mixed phase would be formed in 
the early stage of ultra-relativistic heavy-ion collisions [Q. We have examined the 
surface tension in the flat surface case, as the large flux tube limit, for the various 
Ginzburg-Landau parameter k = m^/\/2mB near the boarder of the type-I and type- 
II dual superconductor. We have found —310 MeV < < 180 MeV for 0.61 < k < 

2.4, when is hxed to = 1.71 GeV. In the type-I dual superconductor with 
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K < l/\/2, the surface tension cToo is positive, so that many flux tubes tend to fuse 
into single huge flux tube [Q. On the other hand, in the type-II case with k > 1/v^, 
(Too is negative, which means that the large size color flux tube is unstable against the 
separation into many thinner flux tubes. The surface tension vanishes, cxoo = 0 just 
at the border between type-I and type-II dual superconductor (k = 1/a/2). 

Finally, we have formulated also the surface tension in a huge color flux tube with 
the large winding number n, which corresponds to the vortex number in the super¬ 
conductivity. We hnd a simple relation between the radius R and the winding number 
n, R = ^/nf, where f is approximately the radius of the unite flux tube. The surface 
tension cr^R) goes down to the asymptotic value cXoo as R ^ oo. Quantitatively, cr^R) 
can be regarded cXoo for R > 2fm. 


H. M. acknowledges the hospitality of the Theory Group at RCNP-Osaka during 
her stay at RCNP, where this study has been worked out. We are grateful to Prof. 
T. Kajino for fruitful discussions on the QCD phase transition in the early Universe. 
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Figure Captions 

Fig.l The effective potential T) as the function of the monopole condensate 

X at various temperatures for case A. The temperature dependence of the 
variable self-coupling A(T) is fixed so as to provide the critical temperature 
at Tc = 0.28 GeV. The absolute minimum at each temperature is denoted 
by the symbol x. At Tc = 0.28 GeV, there appear two absolute minima. 

Fig.2 The effective potential f4fr(x,Tc) at the critical temperature (denoted by 
solid curve) versus the monopole condensate x foi' case A. The dashed curve 
denotes an approximate sine curve. Here, h corresponds to the barrier 
height on the transition between the hadron and the QGP phases, and 
Xh{Tc) denotes the value for the monopole condensate at the non-trivial 
minimum. The energy level is set so as to Ves = 0 at y = 0. 

Fig.3 The surface tension as a function of the Ginzburg-Landau parameter k 

for case B. The surface tension decreases with n. The solid curve is 
obtained by taking the linear temperature dependence for A(T) so as to re¬ 
produce the critical temperature at T = 0.26 GeV, while other parameters 
are fixed. As comparison, the dashed curve is the case with temperature 
dependent v{T), while other parameters as A are fixed to the zero temper¬ 
ature value. 

Fig.4 The profile of the hadron bubble in terms of the monopole condensate 
x{r) near the critical temperature = 0.28 GeV for case A. The solid 
curve denotes the exact solution calculated from T). The dashed 

curve denotes the SG-kink ansatz, 2 xji^ tan“^ where xh is the 

monopole condensate at the non-trivial minimum of Veg{x,T). 

Fig.5 (a) An example of the monopole profile function x(r), approximated by 

the SG-kink curve, as a function of the radial coordinate r. (b) The cor¬ 
responding free energy density. The width of the peak of the free energy 
density corresponds to the thickness of the bubble wall, (c) The total 
free energy of the hadron bubble as the function of various bubble radii 
R. The critical radius Rc is defined at the saddle point of the free energy 
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at each temperature. In this example for case A, is around 26 fm at 
T = 0.2795GeV. 

Fig.6 (a) The surface tension as a function of the critical radius Rc for various 
temperatures, 0.278 GeV < T < 0.28 GeV. The values added nearby the 
curve correspond to the temperatures in unit of GeV. (b) The prohles of 
the monopole condensate at the critical radius Rc at various temperature, 
T = 0.278 GeV, 0.2798 GeV and 0.2799 GeV. The horizontal axis denotes 
the radial coordinate r. The dashed vertical lines show the thickness of 
the bubble wall at each temperature, (c) The solid curve denotes the free 
energy of the hadron bubble with the radius R at various temperatures 
(0.2795 GeV, 0.2797 GeV and 0.2798 GeV) evaluated from the effective 
potential. The dashed curve denotes the free energy approximated by a 
sum of the volume and the surface terms, where ctoo is used. 

Fig.7 The critical radius Rc versus the temperature T for case A. Rc increases 
rapidly with T and goes up to inhnity at the critical temperature Tc- 

Fig.8 The surface tension (Too for the flat interface as a function of the Ginzburg- 

Landau parameter k. This range of k corresponds to 0.5 GeV < ruB < 
2.0 GeV and = 1.71 GeV and \/k^ = 0.42 MeV. At k = 1/v^, the QGD 
vacuum is just at the border between type-I and type-II in terms of the 
dual superconductor. 

Fig.9 (a) The prohle of the monopole condensate x(r) as a function of the radial 
coordinate r for the flux-tube system with various winding numbers n. (b) 
Relation between the flux tube radius R and n (denoted by black dots) 
for K = 0.61 corresponding to the type-I dual superconductor. Here, the 
masses are ms = 2.0 GeV and = 1.71 GeV. This behavior is well 
reproduced by the dotted curve R = f^n, where f = 0.21 fm is almost the 
radius of the single flux tube with n = 1. 

Fig. 10 The surface tension as a function of the radius R of the flux tube. 

The used parameters are the same as in Fig.9. Here, a^o is the value in the 
flat surface case. 
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